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ABSTRACT 

In a recent paper, McLure & Jarvis (2002, astro-ph/0204473, v. 1) reanalyze AGN broad-line reverberation 
measurements presented in Kaspi et al. (2000). They find the broad-line region size and the AGN luminosity 
at 5100 A can be related by Rblr ^ ^5h)0" This differs from the result of Kaspi et al., who found Rblr ^ ^5ioo' 
a departure from expectations for a constant ionization parameter In this note I investigate the source of the 
discrepancy. McLure & Jarvis use quasar luminosities based on the single-epoch, multichannel-spectrophotometer 
measurements of Neugebauer et al. (1987), obtained in 1980. I show that the Neugebauer et al. fluxes are 
systematically higher, by a constant flux offset, compared with the multi-epoch CCD measurements of Kaspi et 
al., obtained concurrently with the echo-mapping radii. In addition, McLure & Jarvis erred when converting these 
fluxes to luminosities. The two effects compound to a typical factor 2 (0.3 dex) overestimate of the luminosities 
of the PG quasars in the sample. Since McLure & Jarvis did adopt the Seyfert luminosities in Kaspi et al. (2000), 
they obtained a slope that is flatter by 0.2. 
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1. INTRODUCTION 

Reverberation mapping of active galactic nuclei (AGNs) ex- 
ploits the light-travel-time delay between continuum variations 
and the response of broad emission lines to measure the sizes 
of AGN broad-line regions (BLRs). This size can then be used, 
together with the estimated velocities of the BLR gas to de- 
rive masses for the central black holes (See Netzer & Peterson 
1996, for a review.) The size-luminosity and mass luminosity 
relations of AGNs may shed new light on understanding these 
objects and the connections between the black holes in AGNs 
and those found in normal nearby galaxies. 

Kaspi et al. (2000) measured reverberation sizes, Rblr, for 
17 PG quasars between 1991-1998 using the Wise Observa- 
tory Im telescope and the Steward Observatory 2.3m telescope, 
with a combination of CCD spectrophotometry and photometry. 
The time averaged flux for each quasar was derived from 20- 
70 epochs per quasar. The rest-frame 5100 A continuum flux 
densities were measured from the ~ 10 A resolution spectra 
while minding potential systematics such as broad-emission- 
line wings and atmospheric absorptions in the higher-redshift 
objects. Following Galactic extinction corrections (small for 
most of the objects) the absolute fluxes were used to calculate 
5100 A luminosities, ALa(5100). The sizes and luminosities of 
Seyfert galaxies with reverberation measurements were com- 
piled from previously published works. A regression analy- 
sis taking into account the errors in both Rblr and ALa(5100) 
showed that Rblr oc [ALa(5100)]° ™±'' °1 This was a surprising 
result, since it has long been speculated that the dependence 
would be to the power 0.5. Such a scaling would lead to an ion- 
ization parameter (ratio of ionizing photon density to electron 
density) at the surface of the BLR clouds that is independent 
of luminosity, and would explain the similarity of AGN spectra 
over many orders of magnitude in luminosity. 

In a recent paper, McLure & Jarvis (2002) examine to what 
degree UV observables, namely 3000 A luminosities and Mg II 
line velocities, can serve the purpose of the optical observ- 
ables - 5100 A luminosities and H/3 widths - used to date in 
AGN black hole estimates. In the course of their work, they re- 
analyze the data presented by Kaspi et al. (2000) and conclude 



that actually Rblr oc [ALa(5 100)]°^"='=° °^ is the best fit, contrary 
to the result of Kaspi et al., and consistent with the expectations 
for a constant ionization parameter. 

Vestergaard (2002) has carried out an analysis along similar 
lines. After studying the various systematics that can affect the 
slope determination, she concludes that the best estimate (her 
equation A5) is Rblr oc [ALa(5100)]°-^^^"°'', consistent with 
the result of Kaspi et al. (2000). The statement by McLure & 
Jarvis (2002), that Vestergaard (2002) found a slope consistent 
with 0.5 but decided to adopt a slope of 0.7 anyway, is incorrect. 

In this Note, I investigate the source of the discrepancy be- 
tween Kaspi et al. and McLure & Jarvis. I show that it arises, 
first, due to the adoption by McLure & Jarvis of old, single- 
epoch, and systematically higher fluxes, but only for the high 
luminosity part of the sample; and second, due to incorrect con- 
version from flux to luminosity. This is confirmed and corrected 
in a revised version of their paper (R. McLure, private commu- 
nication). 

2. ANALYSIS 

McLure & Jarvis chose to base their luminosities for the PG 
quasars in the Kaspi et al. (2000) sample on the measurements 
by Neugebauer et al. (1987). These measurements were ob- 
tained in 1980 with a multichannel spectrophotometer mounted 
on the Palomar 5m telescope. The measurements had coarse 
resolution of ~ 300 A, making the avoidance of emission and 
absorption lines in the measurement more difficult. In the num- 
bers tabulated by Neugebauer et al., interpolation is needed, 
for most of the quasars, to obtain the flux at rest wavelength 
5 1 00 A. As opposed to the measurements of Kaspi et al. (2000), 
only one epoch per object exists, inducing scatter due to the 
time variability of quasars. More importantly, the Neugebauer 
et al. (1987) measurements were obtained about 20 years be- 
fore the Kaspi et al. (2000) size measurements, and it has 
been shown in Seyfert galaxies that BLR size and emission- 
line width change with time in individual objects (Peterson & 
Wandel 2000). 

Figure 1 compares the observer-frame quasar flux densi- 
ties measured at wavelengths corresponding to rest-wavelength 
5100 A, as reported by Kaspi et al. (2000), to the correspond- 
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ing numbers obtained from Neugebauer et al. (1987), after 
the necessary interpolation and conversion to f\ in the latter 
dataset. Apart from the large differences in a few objects, pre- 
sumably due to variability, there is a clear additive systematic 
offset between the datasets, in the sense that the Neugebauer 
et al. (1987) measurements have higher flux. While one can- 
not say for sure which dataset is at fault, it is possible that the 
Palomar data suffer from imperfect sky subtraction - a task that 
is much easier with present-day two-dimesional CCD detectors. 
To get a "third opinion" 1 have compared the Kaspi et al. (2000) 
flux densities to those indicated by the data of Boroson & Green 
(1992). The Boroson & Green (1992) fluxes are generaUy lower 
than the Kaspi et al. fluxes by about 30%. This difference is rea- 
sonable, given that they were obtained through a narrow 1.5" 
sht, and at non-paraUactic angles, and hence some light loss is 
expected. 




fx(5100) [Neugebauer] 



Fig. 1- Comparison of /a [5 100(1 flux densities for PG 
quasars measured by Kaspi et al. (2000) and by Neugebauer 
et al. (1987). Units are 10""' erg s"' cm'^ A"'. 

Fig. 2- Comparison of luminosities \L\ at rest wavelength 
5100 A, calculated from each of the datasets in Fig. 1. 
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Figure 2 examines the effect of this flux discrepancy on the 
calculated luminosities. It compares the 5100 A luminosities 
I have calculated from the Neugebauer et al. (1987) fluxes to 
the luminosities found by Kaspi et al. (2000). In both cases, 
the cosmology assumed is Hq = 75, 0,„ = 1, and iljy = 0. As 
expected, the lower luminosity quasars are most affected, and 
their luminosities are about 0.1 dex higher when based on the 
Neugebauer et al. (1987) numbers. This conclusion was also 
reached by Vestergaard (2002), and Figure 2 is nearly identi- 
cal to her Figure 2a. Indeed, comparing directly her and my 
numbers for the luminosities calculated from the Neugebauer 
et al. data, I find only small differences of order 5%, which are 
ascribable to the interpolation of the Neugebauer et al. data. 

In Figure 3, 1 compare the quasar luminosities listed by 
McLure & Jarvis, to my own calculation based on the same 
Neugebauer et al. data, assuming the cosmology cited by 
McLure & Jarvis, Hq = 70, Qm = 0.3, and Qa = 0.7. Specifi- 
cally, I calculate 

La(5100) = 47rdl(l+z)fx[5l00il +z)], 

where, for the flat (iltotai =1) cosmologies we consider here, 
the luminosity distance is 

Jo 

Figure 3 indicates that there is an error in McLure & Jarvis's 
conversion to luminosities, resulting in too-high luminosities by 
about 0.2 dex, especially at the higher luminosities. R. McLure 
(private communication) confirms that this is a programming 
error that resulted in an extra factor (1 +z)^. 




log(AL^(5100)/Watts) [Mclure & Jarvis] 



Fig. 3- Comparison of 5100 A luminosities reported by 
McLure & Jarvis (2002) with luminosities calculated here 
from the same Neugebauer et al. (1987) data they used, and 
for the same cosmology. 
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Finally, in Figure 4, 1 compare the quasar luminosities listed 
by McLure & Jarvis, to the luminosities of Kaspi et al. (2000), 
after converting the latter to the cosmology and units used 
by the former. Here one sees the compound effect of the 
higher fluxes and the incorrect conversion to luminosity used 
by McLure & Jarvis. This leads to a total offset of about 0.3 
dex (a factor 2). 




log(AL^(5100)/Watts) [Mclure & Jarvis] 



Fig. 4- Comparison of the Kaspi et al. (2000) luminosities 
and the McLure & Jarvis (2002) luminosities, after converting 
the former to the cosmological model and units used by the 
latter. 



The systematically higher, by 0.3 dex, quasar luminosities 
used by McLure & Jarvis, can explain the flatter slope they 
obtained in the Rbim vs La(5100) relation. While preferring 
the Neugebauer et al. measurements over those of Kaspi et al. 
for the quasars, McLure & Jarvis adopted the Seyfert luminosi- 
ties as given in Kaspi et al. (2000). Since the Seyfert Rblr 
and La (5 100) values are, respectively, about 0.8 dex and 1.15 
dex lower than those of the quasars, increasing the quasar lumi- 
nosities by 0.3 dex, will lower the slope of the relation from 
0.7 to 0.55. Some additional flattening also arises from the 
use of the ^A-dominated cosmology, which makes the higher- 
redshift quasars in the sample slightly more luminous. In a re- 
vised version of their paper (R. McLure private communica- 
tion), McLure & Jarvis find a best-fit slope of 0.6 when using 
the Neugebauer et al.(1987) data, and 0.65 ±0.11 when using 
the Kaspi et al. fluxes, fully consistent with Kaspi et al. (2000) 
and Vestergaard (2002). 

I thank Shai Kaspi for stimulating discussions and for sug- 
gesting and performing the comparison to Boroson & Green 
(1992). Ross McLure is thanked for helping sort out the prob- 
lem. 
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